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We show that the periodic orbit sums for 2-dimensional billiards satisfy 
'■Q ' an infinity of exact sum rules. We test such sum rules and demonstrate 

that they can be used to accelerate the convergence of cycle expansions for 
averages such as Lyapunov exponents. 
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1 Introduction 



Periodic orbit theory is a powerful tool for description of chaotic dynamical sys- 
tems ^, pi. However, as one is dealing with infinities of cycles, the formal 
theory is not meaningful unless supplemented by a theory of convergence of cy- 
cle expansions. For nice hyperbolic systems, the theory is well developed, and 
shows that exponentially many cycles suffice to estimate chaotic averages with 
super-exponential accuracy |||, ||, |(|. However, for generic dynamical systems 
with infinitely specified grammars and/or non- hyperbolic phase space regions, 
the convergence of the dynamical zeta functions and spectral determinants cy- 
cle expansions is less remarkable. The infinite symbolic dynamics problem is 
generic, and a variety of strategies for dealing with it have been proposed: sta- 
bility truncations [j?], [|, approximate partitions P], noise regularization |l0| and 



even abandoning the periodic theory altogether [11] 



Computation of periodic orbits for a given system is often a considerable 
investment, as locating exhaustively the periodic orbits of increasing length for 
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flows in higher dimensions can be a demanding chore. It is therefore essential 
that the information obtained be used in an as effective way as possible. Here 
we propose a new, hybrid approach of combining cycle expansions with exact 
results for "nearby" averages, based on the observation that the periodic orbit 
sums sometimes satisfy exact sum rules. 



Studies of convergence of cycle expansions, such as comparisons |12|] of trunca- 
tion errors of the dimension and the topological entropy for the Henon attractor, 
indicate strong correlations in truncation errors for different averages. We pro- 
pose to turn these correlations in our favour, by using the error known exactly 
by a sum rule to improve the estimate for a nearby average for which no exact 
result exist. Billiards provide a convenient, physically motivated testing ground 



for this idea. The approach is inspired by the formula (16) for mean free flight 



time in billiards, so well known to the Russian school that it went unpublished 



for decades [13]. In this paper we show that billiards obey an infinity of exact 
periodic orbit sum rules, and indicate how such rules might be used to accelerate 
convergence of cycle expansions. 

The paper is organised as follows: sect. ^| is a brief summary of the theory 
of periodic orbit averaging. In sect. ^ we review the known exact sum rules for 
billiards, and then generalise them to an infinity of sum rules. In sect. |] we 
present the conventional cycle expansion numerical results for our test system, 
the overlapping three-disk billiard. This system is hyperbolic and does not suffer 
from the intermittency effects that plague billiards such as the stadium and the 
Sinai billiards, but is still "generic" in the sense that its symbolic dynamics is 
arbitrarily complicated. In sect. || and appendix |A| we develop a method which 
utilizes the flow conservation sum rule to accelerate the convergence of cycle 
expansions, and apply the method to our test system. 



2 Periodic orbit averaging 

We start with a summary of the basic formulas of the periodic orbit theory - for 
details the reader can consult refs. |], 

A flow x — > f 1 (x) , x £ Ai , is a continuous mapping f * : Ai — > Ai of the phase 
space Ai onto itself, parameterised by time t. On a suitably defined Poincare 
surface of section V, the dynamics is reduced to a return map 

x^f n {x) (l) 



where n is the "topological time", the number of times the trajectory returns to 
the surface of section. 
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A dynamical zeta function [fl associated with the flow f*(x) is defined as the 
product over all prime cycles p 

l/((z, s, P) = IJ (1 - tp) , * P = tp(*, s, /?) = -L e ^-^^ , (2) 
P I a pI 

where T p , n p and A p are the period, topological length and stability of prime 
cycle p, A p is the integrated observable a(x) evaluated on a single traversal of 
cycle p 



a (f T (xq)) dr (flows) 
A p = { n p -l x e p , (3) 

Y a {f k ( x o)) ( ma P s ) 

fc=0 

where s is a variable dual to the time t, z is dual to the discrete "topological" 
time n, and ip(z, s, /3) is the weight of the cycle p. 

Classical averages over chaotic systems are given by cycle expansions con- 
structed from derivatives of dynamical zeta functions. By expanding the product 
(0) a dynamical zeta function can be represented as a cycle expansion 



t n — t-xi^Zj S, /5) — ( 1) ^Pl^P2 ' ' ' ^Pk 

I 

where the prime on the sum indicates that the sum is restricted to pseudocy- 
cles, built from all distinct products of non-repeating prime cycles weights. The 
pseudocycle topological length, period, integrated observable, and stability are 

A w = A pi + . . . + A Pk , A,,- = A pi Ap 2 • • • Ap k , (5) 

where k = k n is the number of involved prime cycles. For economy of notation 
we shall usually omit the explicit dependence of l/£ and t p on (z, s, (3) whenever 
the dependence is clear from the context. 

Truncation of the dynamical zeta function with respect to the topological 
length n-n < N will be indicated by a subscript 

l/C N (z,s,P) = l- Y! **■ (6) 
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If the system bounded (such that no trajectories escape), the dynamical zeta 
function (|2|) has a leanding zero at 1/(^(1,0,0) = 0. Expressing this condition in 
terms of the cycle expansion (g]) we find that any bound system satisfies the flow 
conservation sum rule (3|: 

l/C(l,0,0) = l-£'(-i)^ +1 -^ = 0. (7) 



The cycle expansions for the phase space average of observable a(x) are given 
by either the integral over the natural measure, or by the cycle expansions 

r ( A )c 

flows: (a) flow = / a{x)p(x)dx = — — (8) 
J M (J- 

f ( A )c 

maps: (a) map = / a(x)p(x)dx = j—^ , (9) 



where p denotes the natural measure. As we shall show in (16) below, the averages 
computed from the two representations of dynamics are related by the mean free 
flight time. 

The cycle expansions required for the evaluation of periodic orbit averages 
(||) and (^) are given by derivatives of the dynamical zeta function with respect 
to (3, s and z 

(A) c = ^l/C(l,0,0)=^'(-l) fc ^ 1 A 7r /|A 7r | (10) 

(T) c = -^i/C(i,o,o) = £'(-i)^/|A w | (ii) 



(n) c = |-l/C(l,0,0) = ^ / (-l) fe '+ 1 n 7r /|A 7r |. (12) 



dz 



3 Periodic obit sum rules for billiards 



We start by reviewing the mean free flight time sum rule for billiards discussed 



by Chernov in ref. [13]. 



In a ci-dimensional billiard, a point particle moves freely inside a domain Q, 
scattering elastically off its boundary dQ. The billiard flow /* on M. = Q x S^ 1 
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(where S is the unit sphere of velocity vectors) has a natural Poincare surface of 
section associated with the boundary 



V = dM = {{q, v) G M : q G dQ and v • n(q) > 0} 



(13) 



where n(q) is the inward normal vector to the boundary at q, defined every- 
where except at the singular set dM* of nondifferentiable points of the bound- 
nary (corners,cusps,etc). In what follows we shall restrict the discussion to (two- 
dimensional) billiards. 

Assume that the particle has unit mass and moves with unit velocity, p\ +P2 = 
1. The cartesian coordinates and their conjugate momenta for the full phase space 
M of the billiard are 

x = (qi,q2,Pi,P2) = (qi,q 2 ,sm(p,cos(j)) . 

Let the Poincare map be the boundary-boundary map / : dM — > dM, and 
parametrise the boundary dM by the Birkhoff (area preserving) coordinates 

x = (s,p s ) , p s = sin 9, 

where s is the arclength measured along the boundary, 9 is the scattering angle 
measured from the outgoing normal, and p s is the component of the momentum 
parallel to the boundary. Both the area of the billiard A = \Q\ and its perimeter 
length L = \dQ\ are assumed finite. Let t(x) be the time of flight until the next 
collision. The continuous trajectory is parametrised by the Birkhoff coordinates 
together with a time coordinate < t < t(x) measured along the ray emanating 
from the boundary point x = (s,p s ). 

The period of a cycle p is the sum of the finite free-flight segments 



where xo = (so,p s ,o) is any of the collision points in cycle p. The mean free 
flight time is the average time of flight between successive bounces off the billiard 
boundary. It can be expressed either as a time average 



n p — 1 



k=0 




i=0 
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or, as a phase space average 



(r) = / r(x)d/i(x), (14) 



V 

where d\i(x) = ds dp s / J-pds dp s is the natural measure. For Hamiltonian flows 
like the billiard flow considered here this is simply the Lebesgue measure. If the 
billiard is ergodic, the time average is defined and independent of xq for almost all 
xq. In order to find an exact expression for the phase space average (r), compute 
the integral over the entire phase space of the billiard, 

/ 5(1 — p\ — pl)dqidq 2 dpidp2 = 2-kA 
JM 

and recompute the same thing in the Birkhoff coordinates, 

/ 5(1 — p\ — p2)dqxdq%dp\dpi = [ dsdp s [ dt = [ r(x)dsdp s 
J Jv Jt=o JV 

= <r> / dsdp s = 2L <r> , (15) 



where L is the circumference of the billiard. Hence the mean free flight time is a 
purely geometric property of the billiard, 

(t) = IT j, (16) 

the ratio of its perimeter to its area. The relation is a consequence of the Liouville 
measure being constant and apply to any billiard regardless of whether its phase 
space is mixed or not. For ergodic systems the periodic orbit theory gives a cycle 
expansion formula @ for the mean free flight time 



(r) = (17) 



If we know (r) this formula enables us to relate any discrete time average @ 
computed from the map to the continuous time averages @ computed on the 
flow. They are linked by the mean free flight time formula 

( a ) map = ( a ) flow ( T ) ■ ( 18 ) 
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In what follows we will restrict our attention to map averages, and omit the 
subscript {...) map -»(...>• 

As the next example of a periodic orbit sum rule, consider the case of the 
observable being the transverse momentum change at collision, a = 2cos#. The 
corresponding sum rule is called the pressure sum rule because it is related to the 
pressure exerted by the particle on the billiard boundary. 

The average pressure is given by the relation P = F/\dQ\, where F is the 
time average of momentum change, that is the force the particle exerts against 
the boundary. The momentum change per bounce equals twice the transverse 
momentum at the collision, so the average force per bounce is 

(F) m a P = I 2p±(xW(x) = f J, / f ^ 2cos0 casOdOds = ~ . (19) 
y Jv Jp dsdp s Jqq 7-71-/2 2 

Hence the pressure for a flow becomes 



The exact averages (16), (EH) apply to billiards of any shape, ergodic or not. 
As both the mean free path and pressure can be calculated by means of cycle 
expansions, these relations leads to exact billiard sum rules. 

Now we note that as the Liouville invariant measure of the map is constant, 
any average of an observable a(x), defined in terms of dA4 coordinates x = (s,p s ) 
can be expressed in terms of a simple integral. For each such observable we obtain 
an exact periodic orbit sum rule 

(a) = = / a(x)d^x). (21) 
(n) ( Jv 

Surprisingly enough, this uncountable infinity of sum rules seems not to have 
been noted in the literature. 

The formula (21) does not allow for analytical computation of every average 
we want to compute in a billiard. Consider the simplest nontrivial average worth 
study in billiards, the diffusion coefficient 



i i d 2 

2d (n} c df3 2 



(^) = ^7ZV-7^( 1 /C)l /3=0 - (22) 
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It requires evaluation of a second derivative of the relevant dynamical zeta func- 
tion, this means that two-point correlations of the observable along cycles will 
enter the averaging formulas, so the average can not be computed from one iterate 
of the map. 

Another quantity of interest is the Lyapunov exponent. Let A(xo,n) is the 
largest eigenvalue of the Jacobian of the nth iterate of the map. The (largest) 
Lyapunov exponent is defined as 

A= lim —log\A(xQ,n)\. 

n— »oo ji 

The cycle expansion formulas in sect. [2] compute 

A= lim - / dp( x) log \A(x,n)\, (23) 

n-*oo n J 

that is, a combination of time and phase space averages. Note that if A{xo,n) is 
multiplicative, A(xo,n) = ]lfc=o ^~(f k ( x o))> then the integral in (|23|) is indepen- 
dent of n: setting n = 1, reduces the problem to one iterate of the map. This is 
the case for 1-dimensional maps. However in most cases the invariant measure 
p(x) is not known a priori, and there is no simple exact formula for the average. 

For billiards the invariant density is known but the expanding stability eigen- 
value A(xo,n) is not multiplicative along an arbitrary trajectory, , and the integral 
in (^) is dependent on n. It is possible to derive a multiplicative evolution op- 
erator for this purpose [fbf l. However, for the purpose at hand the naive cycle 
expansion formulas still apply, because A(xQ,n) is multiplicative for repeats of 
periodic orbits. By defining the cycle weight 

= \A p f 



the cycle expansion for the Lyapunov exponent is given by 
(ln|A|), 

<A> = (24) 

So even though Lyapunov only requires computation of two first order deriva- 
tives of the dynamical zeta function, it requires n-point correlations to all orders 
and cannot be computed by a sum rule. 



In the case when ( |18| ) the Lyapunov exponent of the flow to the Lyapunov 
exponent of the corresponding Poincare return map, the relation is known as the 



Abramov's formula [15| 



S 



Figure 1: The overlapping 3-disk billiard. A point-like particle moves inside the 
billiard bouncing specularly off the boundary. Shown is a cycle of topological 
length 4. 



4 Three overlapping disks billiard 



We will test the above sum rules on cycle expansions for a concrete system, 
the overlapping 3-disk billiard. This billiard consists of three disks of radius 
a centered on the corners of an equilateral triangle with sides R. There is a 
finite enclosure (see fig. |]) between the disks if v3 < R/a < 2. This enclosure 
defines the billiard domain Q. One of the limits R/a — > y/3 corresponds to 
the integrable equilateral triangle billiard. The other limiting case R/a = 2 
exhibits intermittency with infinite sequences of periodic orbits whose periods T p 
accumulate to finite limits, and where stabilities fall off as some power , where 
n p is the topological length. 

The C^ v symmetry of the billiard enables us to work in the fundamental 



domain [16]. The fundamental domain is one 6th slice of the billiard domain, 
fenced in by the symmetry axes of the billiard. In what follows we are only 
interested in the lowest eigenvalue and therefore we restrict our computations to 
the fully symmetric A\ subspace. The fundamental domain symbolic dynamics 
is binary, but is not of the finite subshift type; its full specification would require 
an infinity of pruning rules of arbitrary length. 

The mean free flight time (|l6|) for the overlapping 3-disk billiard can be found 
by geometric considerations, 



^^ vw-M-w (25) 



where r = y/a 2 — (R/2) 2 and 9 = 7r/6 — arcsin(r/a). We shall set a = 1 through- 
out this paper, and parameterise the billiard by the center-to-center distance R. 
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R 


(r) (A) Anum L (P) 


# cycles 


1.85 
1.9 


0.102 0.523 1.57 
0.1401 0.6036 0.60363 1.570 


342 
525 



Table 1: The mean free flight time (t) , the average pressure P , and the best estimate of 
the Lyapunov exponent A computed by cycle expansion as function of the 3-disk center- 
to-center separation R used in our numerical tests, with disk radius fixed to a = 1. For 
R = 1.9 a numerically computed Lyapunov reference value obtained by direct simulation 
using 10 10 bounces is displayed. The total numbers of the fundamental domain prime 
cycles used in the cycle expansion computations are also indicated. 




2 4 6 8 10 12 
N 



Figure 2: Convergence of cycle expansions: Deviation of cycle expansions trun- 
cated to the topological length N from exact sum rules for (o) flow conservation 
([?]) and (□) mean free flight time ([T^), 

All our numerical tests are done for R = 1.9. Results for this parameter value, 
as well as for R = 1.85 are shown in table |l[ 

Fig. [2] illustrates the convergence of finite topological length cycle expansions 
for the flow conservation sum rule (0) and for the mean free flight time sum rule 
(p^). As the exact result is known, we plot the logarithm of the error as function 
of the truncation length N. 

The overall exponential convergence reflects the existence of a gap, the dyn- 
amical zeta functions are analytic beyond z = 1. The "irregular" oscillations in 
fig. ^ are typical for systems with complicated symbolic dynamics. For systems 
with finite subshift symbolic dynamics the oscillations ceases when the cutoff N 
exceeds the longest forbidden substring and if the full spectral determinant is 
used, super-exponential convergence sets in 

One should note the coincidence of the peaks and dips of the two curves. 
This type of correlations between coefficients between different power series will 
be important in the following. 
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5 Utilizing exact sum rules 



We shall illustrate the utility of exact sum rules in accelerating the convergence 
of cycle expansions by applying flow conservation sum rule to the problem of 
computing the mean free flight time fll6|). As we already have the exact formula 
for this average, we will be able to compute the exact error in the various estimates 
and compare them. We will then apply the same technique to evaluate the 
Lyapunov exponent, for which no exact formula exists. 

The idea is to use the flow conservation sum rule (0) to improve the numerator 
and denominator of (|l7|) separately. 

We begin by the denominator. The general problem is to find a good estimate 
of the derivative F'(zq) at the first root F(zq) = of a function given by a power 
series F(z) = Y^k=Q^k zk where the coefficients are known only up to order N 



N 

F {N) (z) = Y / hz k , (26) 

k=0 

and we know (appendix |A|) that a good estimate is given by 

TV + 1 

F'(zo) « F N) (z ) - —^F (N) (z ). (27) 

We argue in appendix |A| that the error in the above estimate is suppressed com- 
pared to the error of the estimate F! n Jzq) by a factor q whose asymptotic be- 
haviour is 



?~1/JV. (28) 
The improved estimate of the derivative of the zeta function ( [12] ) thus reads 
(rc>C,acc « (n)c,(N) - ( N + 1 )C{n ) (1,0,0). (29) 

The continuous time average by invoking the sum rule ([?]) is similar to the 
previous one but l/£(l,s,0) = is now a Dirichlet series in s. The basic idea is 
to start by expanding the zeta function around some point s = sq 



kjr 1 



l/C(M,0) = l-E / ( H | e 
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= i-j2' t*)—t e -«T. g (*o-*) fc ^ 

Since we have only a finite number of of pseudocycles at our disposal, it is not 
meaningful to consider high powers k in the above Taylor series. So we truncate 
the series 

VWM.0) - - E' tgf^E**^ 

7r:n 7r <AT ^""l k=0 



l-E^o-f. (30) 



fc=0 



where 

T k\ = y> ( _i)^-i e is 

is the fcth moment of the average cycle time. The sets of periodic orbits going 
into the calculation are still being truncated according to their topological length. 
This finite Taylor series is the analogue of the truncated function (z) treated 
above; s = corresponds to z = 1, and s = so to z = 0. We can use (p7| ) and 
write down the improved estimate 

(T) c , acc « (T) C)(JV) - ^±C^ >M) (1,0,0) (31) 

However there is now an additional complication due to the fact that not all 
available (T k ) are exact. So how should we choose sq and Ml We choose sq to 



C 

lie somewhere in the range 1/T max < sq < 1/T m j n where the 

Tmin and T max are 

the smallest and largest period in the sample for a particular topological length 
cutoff N. 

The next question is, for a given sq, what is the number of accurate coefficients 
(r k ^ ^ /k\l We see from ( ^0[ ) that pseudocycles are suppressed with their 

length according to the function T k exp(— sqT) having its maximum at T = k/sQ. 
So the coefficients with k soT max can be expected to be accurate. However, 
as the majority of cycles have periods close to T max we want to make use of the 
information they carry. We have found it preferable to include a large number of 
fairly accurate coefficients rather than a small number of very accurate ones. So 
we choose the maximum power M to be given by the average cycle length 

M = s T p 

p n„=N-l 
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1 1.5 2 

InlNI 



Figure 3: The error suppression factor for (□) maps and for (o) flows, applied to 
the 3-disk system with R = 1.9. Here we have used an extrapolated value from the 
cycle expansion as the best asymptotic estimate. Both error suppressions display 
the estimated error decrease and demonstrate that the sum rules do improve 
convergence. 



The error of the improved estimate is suppressed compared to the error of 
the traditional estimate by factor we call q, see appendix |A]. This (/-factor is 
plotted in fig. ||. It decreases (apart from oscillations) as the estimated iV -1 
error suppression derived for maps. 

The calculation of the integrated observable amounts to evaluating the (3 
derivatives of the dynamical zeta functions. The role of [3 is completely analogous 
to that of s. With (3 viewed as a complex variable, the dynamical zeta function 
1/C(1,0, (3) is a Dirichlet series in (3 and the above methods can be used to 
compute ^1/C(1,0, 0). Here similar criteria apply to (3q and N as for (^Tj) : f3o 

close to 1/A m i n and M = (3q A p 



5.1 Improvement on the averages 

So far we have improved the numerator and the denominator of ([l?]) and (||) 
separately. The errors of both are suppressed by a factor q 0(1/N) compared 
to unaccelerated estimates. We have also seen (fig. |3|) that, both before and after 
resummation, their behavior versus the cutoff N are highly correlated. So it is 
not obvious how the resulting average should be improved, indeed it is not clear 
whether it is improved at all. 

The accelerated cycle expansion for an observable a(x) using our method is 



( a )acc = -J-T . ( 32 ) 
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Figure 4: The error suppression factor d33| ) for : (O) The accelerated mean 
free flight timesum rule, (o) pressure sum rule. (□) The Lyapunov exponent. 
At topological length 12 the accuracy of the accelerated Lyapunov exponent has 
reached the best estimate from direct numerical simulation (see Table 1). 



and the error suppression, the g-factor for the observable a{x) is 



( a )acc ( a )exa,l 



(a) - (a) 



(33) 



exact 



We compute this q factor numerically for three different averages: 



(i) The mean free flight time (r) by (fHj). The exact result is given by (16). 

(ii) The average force {F) map by (|l9D , associated with the pressure sum rule. 
The exact result is given by (|20|), 



(iii) The Lyapunov exponent by (|24|). The reference value of the Lyapunov 
exponent is obtained by numerical simulation, see Table 1. 



The results are summarized in fig. ||. The accelerated cycle expansions are 
clearly better than the standard cycle expansions. The error suppression factors 
appear to decrease exponentially, and therefore the acceleration techniques has 
for the 3-disk system increased the correlation between the expansions leading to 
a faster convergence for the averages. 



6 Conclusion 



In this paper we have achieved two objectives: (i) We have derived an infinite 
number of exact periodic orbit sum rules for billiards (pl|). Such sum rules enable 
us to make exact computations of some statistical averages for billiards, such as 
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the mean free flight time (|ilf ) and pressure (|20|), (ii) We have derived the im- 
proved estimate ( |3l| ) which combines the flow conservation sum rule (|7[) with 
the cycle expansions. In order to measure the convergence acceleration, we have 
introduced the error suppression factor ( |33|) that gauges the improvement of the 
accelerated cycle expansions relative to the unaccelerated ones. We thus demon- 
strate that exact sum rules can be used to accelerate convergence for observables 
for which no exact results exist, see fig. [|. 

A challenge for the future is to utilize such infinities of sum rules for billiards 
in the classical applications (other than the Lyapunov exponent studied here), as 
well as in the semi-classical applications of periodic orbit theory. 

This work was supported by the Swedish Natural Science Research Council 
(NFR) under contract no. F-AA/FU 06420-312 and no. F-AA/FU 06420-313. 
PD thanks NORDITA for partial support. SFN thanks PD and KTH for hospi- 
tality. 



A Resummation of power series 



Consider a function F(z) given by a power series, where only a finite number of 
coefficients are known. 

N 

F {N) (z) = J2hz k . (34) 

k=0 

We assume that F(zq) = for some zq and we wish to estimate the first 
derivative F'(zq) (and possibly higher derivatives) as accurately as possible. The 
general problem is to transform the Taylor series around z = into a Taylor 
series around z = zq, and extract the desired coeficients. This is done by the 
ansatz 

N N+l 

£ b k z k = J2 <h(z - z y + 0(z N+1 ) (35) 

k=0 1=1 

Note that the sum rule is built into this ansatz by setting ao = 0. Our aim is to 
determine F'(zq) = a\. We keep the number of known and unknown coefficients 
equal so that the system of equations is solvable. 



Expanding the right hand side ( p5| ) binomially 

(36) 




i=0 t=l j=0 
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we obtain the linear system of equations 



N+l / A 



b i = H ■ {-Vf~ j a h 1 < i < N+ 1 and < j < N. (37) 

max(j,l) W 



In order to transform this to matrix equations where all indices range from 1 to 
n = N + 1 we define vectors 

(b)i = «j -1 6i_i, (a)i = 4«i , (38) 
and rewrite (|37|) as 



b = Ma, (M) 4J = .M, l<i,i<n(-iy- m . (39) 



We use a convention that (",) = if m is out of range. This system may readily 
be solved. Define the matrix L by 



W« = \ n ■ ^ (40) 



~ \ i < j 

Then 



(LM)jj = (-l)W 1 P _ M _> (LM): 1 = - ( 3 _ \ j , (41) 



and the explicit solution is 

a=(LM) x Lb (42) 
In particular 



(a)i = -n(b)i-(n-l) (b) 2 - . . . - 1 (b) n = - - A; + l)(b)* 

k=l 

N 

z oai = (a) 1 = -Y f (n-k)4h = z Fl N) (z )-(N + l)F w (zo) (43) 

fc=0 
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So our improved estimate of F'(zq) is 

F'(z ) « ai = F^) (zq) ~(N + l^ify) (z ) (44) 

The error is suppressed by a factor 

= F'^-nz, 1 F N (z )-F'(z ) = ^Jk - n)b k 
q F' N (z )-F'(z ) E%=nkh 1 ' 

To get this on a more handy form we use summation by parts, that is, we define 

oo 

l/g=l+ v oT" - , s k = Y,bj (46) 

If F(z) is the spectral determinant for a <i-dimensional Axiom-A map the 
coefficients of the power series expansion are super-exponentially bounded 



C a A- a nW/d < \b n \ < C b A^ 1+1/d (47) 



where 1 < A& < A a . Assuming moreover that the signs of the coefficients settle 
down to some periodic pattern, one can show that the error supression factor has 
the following asymptotic behaviour 

q ~ n"( 1+1 ^ (48) 



In this paper we focus on a hyperbolic systems whose symbolic dynamics can- 
not be finitely specified. In that case the bound on the coefficients is exponential 

C a A~ n < \b n \ < C 6 A 6 -", (49) 



and nothing can be said about the signs, as they can oscillate in a completely 
irregular fashion [17]. It seems difficult to obtain proper bounds on q in a general 



setting. In the case at hand we can only provide a qualified guess on the decrease 
of the error supression factor 

q ~ n- 1 (50) 
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Figure 5: Error suppression factor ( |33| ) (y-axis) versus truncation in topological 
length N(x-axis) for the tent map (with a "typical" slope value A). 



Some evidence for this behavior can be provided by the tent map 



F(x) 



Ax for x < \ 
-A(x — 1) for x > | 



(51) 



The expansion rate is uniform but complete symbolic dynamics is lacking in the 
generic case. In fig. || we plot the q-factor for the tent map for a randomly chosen 
parameter versus N. It conforms with the predicted 1/N behavior. 

The ansatz (^) used here is the simplest conceivable and it led to very simple 
formulas. The only requirement is that the dynamical zeta function is analytic 
in a disk z < R where R > 1. This excludes strongly intermittent systems where 
a more refined ansatz is needed ]l8|. If one has some explicit knowledge of the 
nature of the leading singularity of the dynamical zeta function, one can taylor 
a more specific ansatz. 
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